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Abstract 

Pseudoscalar particles (j) usually couple electromagnetically by an interaction of the form 
\g(t)FF, allowing them to convert to photons in the presence of magnetic fields. Notably, new 
low-mass pseudoscalars emitted from supernova (SN) 1987A would have been converted to 7-rays 
in the intervening magnetic field of the galaxy. Therefore, measurements by the Solar Maxi- 
mum Mission (SMM) Gamma-Ray Spectrometer (GRS) can limit the inverse coupling constant 
to g-^ > IxlO" GeV, assuming the pseudoscalar is massless. This is an improvement over other 
astrophysical limits of a factor of about 2.5. 
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I. INTRODUCTION 



Light neutral pseudoscalars arise naturally as a result of spontaneously broken global symme- 
tries ||l|. They will be truly massless if the symmetry is not anomalous. Such particles generically 
couple to the electromagnetic field through an interaction of the form 

L = IgcPF^^F^, , (1) 

where F^i, is the dual of the electromagnetic field-strength tensor F^'^ and g is a coupling constant 
of dimension (energy) which is related to the scale of symmetry breaking. This coupling leads to 
the interconversion of photons and pseudoscalars in external electric or magnetic fields. Just as the 
pion decay constant was originally measured by this Primakoff effect, it is natural to investigate 
the existence of new pseudoscalars by this method. For example, two beautiful experiments to 
search for the presence of dark-matter axions in our galaxy are currently under way ||2|. 

If new low-mass bosons exist, stars will be powerful sources for their production |3|,^. There- 
fore, a classic way to constrain the coupling g is from observational limits on anomalous stel- 
lar energy losses. The observed helium-burning lifetime of horizontal-branch stars yields a limit 
> l.TxlO^*^ GeV, applicable if the particle mass is not much larger than the stellar core tem- 
perature of about lOkeV and thus covers the important case of axions 0j5[. 

However, if the pseudoscalar mass is very small, more restrictive limits can be obtained from 
the possible conversion of stellar particle fiuxes in the large-scale magnetic field of the galaxy. 
This effect would lead to apparent x-ray or 7-ray fiuxes from stars since the emitted pseudoscalars 
have energies representative of the interior temperature. Essentially, then, one would observe 
photons from the stellar core which were converted there to pseudoscalars in the electric fields 
of the charged particles of the heat bath, and which were then back-converted into photons by 
the galactic magnetic field. One of us (E.G.) has previously applied this argument to the red 
supergiant a-Ori (Betelgeuse) and found a limit g^^ > 4x10^° GeV |6|. Later, one of us (G.R.) 
pointed out that this method could also be applied to the pseudoscalar flux from supernova (SN) 
1987A Q. The Gamma-Ray Spectrometer (GRS) on the Solar Maximum Mission (SMM) satellite 
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has set very restrictive limits on a possible 7-ray burst in conjunction with the observed SN 1987A 
neutrino signal, a fact which has been used to derive limits on neutrino radiative decays Q. We 
presently show that these observations yield a constraint g^"^ > 1x10^^ GeV, more restrictive than 
all previous bounds. Of course, it will apply only to massless pseudoscalars (true Goldstone bosons) 
and thus excludes axions for which the horizontal-branch star limit remains the most restrictive. 

In Section 2 we discuss the theory of pseudoscalar production in the supernova's interior. In 
Section 3 we discuss their conversion to photons in the presence of the galaxy's magnetic field 
and the limit which such conversion can achieve. Observational limits on converted photons by the 
Solar Maximum Mission satellite is discussed in Section 4. We summarize our findings in Section 5. 

II. SUPERNOVA PRODUCTION OF PSEUDOSCALARS 

Pseudoscalars with a coupling of the form (|^) are produced primarily by the Primakoff process. 
In vacuum, its cross section is logarithmically infinite due to the infinite range of the Coulomb 
potential. In a medium, however, Debye screening cuts off this divergence and leads to a conversion 
rate per unit time for photons to pseudoscalars of Q 

(2) 

where E is the photon energy, T the temperature, and k the inverse Debye screening length. The 
overall factor of in this equation is simply a convenient way of writing the density of scattering 
targets, and has nothing especially to do with Debye screening. The other appearances of k? arise 
from the finite range of the electric field surrounding charged particles in the plasma. 

In order to derive an expression for n which is appropriate for the conditions of a SN core, 
we note that the relevant charged particles are equal number densities of electrons and protons. 
However, because of their large mass difference the electrons are very degenerate, while the protons 
are nearly nondegenerate, or at most partially degenerate. Therefore, electrons are essentially 
unavailable as scattering targets because their phase space is almost completely Pauli blocked. 
Further, their degeneracy causes them to form a "stiff" background which is difficult to polarize 
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so that they can be ignored for screening as weh. Because of the large temperature, the plasma 
coupling parameter for the protons is much smaller than unity in spite of the large density so that 
the plasma is weakly coupled in the electromagnetic sense. Therefore, k is given by the Debye 
formula 

^2 ^ l^j^^ ^ (3) 

where Up is the number density of the protons and a is the fine-structure constant. 

Multiplying Eq. (^) with the density of thermal photons, we find the pseudoscalar volume 
production rate per unit energy to be 

where = kP/AT'^. We stress that it is possible to treat photons essentially as massless particles 
because the plasma frequency for the conditions at hand is small compared with the temperature. 

In order to calculate the total expected flux we use the numerical SN model S2BH_0 of Keil, 
Janka, and Raffelt which is a representative case. We had numerical details at hand for t = 1, 
5, and 10 seconds after core bounce. The values of T and ^'^ as a function of radius are plotted in 
Fig. 1. Further, in Fig. 2 we show T/Ep, where the nonrelativistic proton Fermi energy is given 
by Ep = pp/2mj\(. (The Fermi momentum is related to the proton density by Hp = pp/Svr^.) 
Degeneracy effects typically become truly important only for Ep/T ^ 3 so that our approximation 
of treating protons as nondegenerate is quite reasonable. 

We then integrate over the volume of the star to get the total rate dNfp / dE at which particles 
are produced. It is graphed in Fig. 3 for a coupling strength of g^^ = 10^'^ GeV. For the earliest 
time, t = \ sec, we may be slightly overestimating the production because of the proton degeneracy. 

Once the pseudoscalars are produced, they will escape the supernova provided their mean free 
path A for backconversion exceeds the size of the SN core. The backconversion rate is just twice 
the rate given in Eq. (^). Typical photon and pseudoscalar energies will be around E = 3T, so 
that typically E/2k, = 3/^. In this range the square bracket in Eq. (^) is of order unity. Therefore, 
A^-*^ ~ g'^K?T/167r, almost independently of the energy. If we use T = 30MeV and k = 50MeV as 



characteristic values, we find A ~ g^^ 10^^ cm where g^Q = g~^ /IQ^^ GeV. Therefore, in the range 
of couphng constants which is of interest here, the pseudoscalars escape freely from the SN core 
once produced. 

III. CONVERSION IN THE GALACTIC MAGNETIC FIELD 

Having produced our pseudoscalars, we now wish to detect them. To do so, we will back- 
convert them in the galactic magnetic field. In a field which is roughly homogeneous on scales much 
larger than the pseudoscalar wavelength, this conversion process can be viewed as an oscillation 
phenomenon much like that of neutrino oscillations [|lO|. However, oscillation phenomena will 
manifest themselves only if the momentum transfer (m^ — m'^)/2E which is needed to convert a 
pseudoscalar into a photon exceeds the inverse scale over which the i?-field is roughly homogeneous. 
For photons, the "mass" is given by the plasma frequency = Airane/me where He of order 
0.03 cm~'^ is the interstellar electron density. It corresponds to = 0.6x10^^^ eV or, with a 
photon energy E = 100 MeV, to my2E = 2x10"^^ eV = (30Mpc)~^ The relevant spatial extent 
of the magnetic field will turn out to be of order 1 kpc so that photons can be viewed as massless 
for the present purposes in spite of the plasma effect. Similarly, for the pseudoscalars to count as 
effectively massless we need to require (m^/2£')^^ > Ikpc or m,^ < 10~^ eV. 

If photons and pseudoscalars are effectively massless in this sense, the conversion probability 
is given by 

P = \g^Bie , (5) 

where -B_l is the magnetic field perpendicular to the line of sight, and i is the distance over which 
the magnetic field is effectively constant. 

The magnetic field of the galaxy has considerable structure. There is certainly a toroidal 
component with a magnitude of about 2 /zG and a coherence length of about 1 kpc. There are also 
other components, which may be roughly characterized as a random contribution of magnitude 
5 /iG with a coherence length of perhaps 10 pc. Because the random component cannot contribute 
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coherently, we anticipate that the toroidal field will be the dominant contribution. Hence we will 
assume that there is a constant field of 2 fxG. Since 1 kpc is also the approximate thickness of the 
disk where we expect this galactic magnetic field, we will assume a constant field on this scale, 
which cuts off suddenly after a distance of 1 kpc. There may also be extragalactic fields between 
us and the Large Magellanic Cloud; if so, they will only increase the effect. Obviously, there is 
considerable uncertainty in this calculation, and hence our final result can not be taken as precise. 

SN 1987A is at a galactic latitude b = —32.1° and longitude I = 279.6°. This unfortunately 
means that we are looking primarily along the direction of the magnetic field, so that 

P= Ig'^B^fil-cos'^bsm'^l) = l0.30 g'^B^f . (6) 

The photon flux at Earth is then given by 

dE dE 4ttD^ ' ^ ^ 

where i = 1 kpc is the length of the conversion region, D = 50 kpc is the distance to SN 1987A, 
and Bl = 0.30 (2 fiGf. For a couphng g-^ = 10^°GeV, this yields the fluxes plotted in Fig. 3. 



IV. DETECTION OF GAMMA-RAYS 

Not surprisingly, no detector was oriented speciflcally in the direction of SN 1987A at the time 
of core collapse. Fortunately, the Gamma Ray Spectrometer (GRS) on the Solar Maximum Mission 
(SMM) satellite, though pointed in the direction of the Sun, was still able to detect gamma rays 



coming from the direction of SN 1987A |11]. The precise timing of the core collapse is known 



from the observation of neutrinos by the Irvine-Michigan-Brookhaven |12| and Kamiokande II |13] 
detectors. Because electron neutrinos and photons are both effectively massless they should arrive 
simultaneously at Earth. Thus we know exactly which period to search for an excess of photons 
from SN 1987A. 

During the 10.24 seconds that the neutrinos were detected from SN 1987A, the GRS searched for 
photons in three energy bins, 4.1-6.4 MeV, 10-25 MeV, and 25-100 MeV. The 95% confidence limits 
on the total fluence of photons during this time period were 0.9, 0.4, and 0.6 7cm~^ respectively 
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|11|. Because of the rising shape of our spectrum, the best hmit for us comes from considering 
the highest of these bins. The detailed hmit depends on the shape of the spectrum, and the hmit 
cited in |11] assumes a spectrum which faUs as E~'^. Without knowing more about the energy 
dependence of the detector's response, we cannot get a precise limit, but we will assume that the 
limit 0.6 7cm^^ applies. 

In the highest energy bin, the 7-ray flux from SN 1987A can be found by integrating the 
rates plotted in Fig. 3. For a coupling of g^^ = 10^'^GeV, this turns out to be 2880, 2120, and 
11807cm~^s~^, at t = 1, 5, and 10 seconds, respectively. Recall that at t = 1 second, there should 
be a noticeable reduction of the rate from that calculated due to Pauli blocking of the final-state 
protons. We expect this blocking to be even larger at earlier times, so we ignore all contributions 
from t < 1 second, and assume that our calculations are accurate over the rest of the time period. 
A linear fit gives a fiuence (time-integrated fiux) of 1.8x 10'* 7 cm~^ for the period from t = 1 to 10 
seconds. 

Nominally, this implies a limit of g^^ > 1.32x10^* GeV. However, there are considerable un- 
certainties involved in the calculation, principally in the details of the galactic magnetic field, and 
the energy dependence of the detector response. If we assume that there is an uncertainty of a 
factor of three in our calculated fiuence, then our limit is weakened to g^^ > 1.0x10^* GeV. Note 
that since the photon fiux is proportional to g'^, the g limit is rather insensitive to small errors in 
the fiux or detection calculation. 



V. DISCUSSION AND SUMMARY 

We have found that the absence of excess 7-rays from SN 1987A limits the photon coupling of 
a nearly massless pseudoscalar g^^ > 1x10*^ GeV. This limit is significantly more restrictive than 
g^^ > 1.7x10^'^ GeV which is found from the classic energy-loss argument applied to horizontal- 
branch stars However, the latter argument is valid for pseudoscalar masses of up to about 
lOkeV and thus covers the important case of axions. Of course, for axions even more restrictive 
limits obtain from their coupling to nucleons on the basis of the energy-loss argument applied to 
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SN 1987A. Our present result is primarily useful for hypothetical massless Goldstone bosons which 
do not exhibit significant couplings to "normal" fermions such as electrons or nucleons. If they 
had such couplings one would typically expect that these couplings yield more restrictive limits on 
the underlying scale of symmetry breaking. Ignoring such couplings also justifies that we focused 
on the Primakoff effect to estimate the mean free path in the SN core. In principle, of course, 
arbitrary pseudoscalars could be trapped by interactions with nucleons or electrons. 

In a previous paper one of us (E.G.) considered what limit had been set or could be set 
by looking for 7-rays from stars, particularly the nearby supergiant a-Ori (Betelgeuse). A limit of 
g^^ > 4x10^^ GeV was set by reexamining old data, and it was pointed out that this limit could 
be improved to the level of 1 x 10^^ GeV if a dedicated observation were made. We have achieved 
this goal with the existing SN 1987A data. 
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Figure Captions 



Fig. 1. The temperature T (la) and the parameter ^ (lb) as functions of radius for times t=l, 5, 
and 10 seconds, based on the model S2BH_0 of Keil, Janka and Raffelt |^. 

Fig. 2. The parameter T/Ep (proton Fermi energy Ep) as a function of radius r for t=l, 5, and 
10 seconds after core bounce, based on the model S2BH_0 of Keil, Janka and Raffelt [^. 

Fig. 3. The number of pseudoscalars produced (left scale) and resulting 7-ray photon flux at the 
Earth (right scale) for a coupling g^^ = 10^'^GeV for times t=l, 5, and 10 seconds. 
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